Evolution of transonicity in an accretion disc 
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Abstract. For inviscid, rotational accretion flows driven by a general pseudo- 
Newtonian potential on to a Schwarzschild black hole, the only possible fixed points 
are saddle points and centre-type points. For the specific choice of the Newtonian 
potential, the flow has only two critical points, of which the outer one is a saddle point 
while the inner one is a centre- type point. A restrictive upper bound is imposed on the 
admissible range of values of the angular momentum of sub-Keplerian flows through a 
saddle point. These flows are very unstable to any deviation from a necessarily precise 
boundary condition. The difficulties against the physical realisability of a solution 
passing through the saddle point have been addressed through a temporal evolution of 
the flow, which gives a non-perturbative mechanism for selecting a transonic solution 
passing through the saddle point. An equation of motion for a real-time perturbation 
about the stationary flows reveals a very close correspondence with the metric of an 
acoustic black hole, which is also an indication of the primacy of transonicity. 
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1. Introduction 

In accretion processes critical (transonic) flows — flows which are regular through a 
critical point — are of great interest [T|. Transonicity would imply that the bulk velocity 
of the flow would be matched by the speed of acoustic propagation in the accreting fluid. 
In this situation a subsonic-to-supersonic transition or vice versa can take place in the 
flow either continuously or discontinuously. In the former case the flow is smooth and 
regular through a critical point for transonicity (more specially, this can be a sonic 
point, at which the velocity of the bulk flow exactly matches the speed of sound), while 
in the latter case, there will arise a shock [2] . The possibility of both kinds of transition 
happening in an accreting system is very much real, and much effort has been made so far 
in studying these phenomena. For accretion on to a black hole especially, the argument 
that the inner boundary condition at the event horizon will lead to the exhibition of 
transonic properties in the flow, has been established well pQ. 

A paradigmatic astrophysical example of a transonic flow is the Bondi solution in 
steady spherically symmetric accretion [HE]. What is striking about the Bondi solution 
is that while the question of its realisability is susceptible to great instabilities arising 
from infinitesimal deviations from an absolutely precise prescription of the boundary 
condition in the steady limit of the hydrodynamic flow, it is easy to lock on to the 
Bondi solution when the temporal evolution of the flow is followed |4J. This dynamic 
and non-perturbative selection mechanism of the transonic solution agrees quite closely 
with Bondi's conjecture that it is the criterion of minimum energy that will make the 
transonic solution the favoured one [21 [5] . 

The appeal of the spherically symmetric flow, however, is limited by its not 
accounting for the fact that in a realistic situation, the infalling matter would be in 
possession of angular momentum — and hence the process of infall should lead to the 
formation of what is known as an accretion disc, i.e. the accretion process would be 
axisymmetric in nature. A substantial body of work over the past many years, has 
argued well the case for transonicity in axisymmetric accretion [H IS, El El EJ [10, [HI 
[T2l [T3| HU [15j [16j [T7J [18]. Since in either of the two kinds of astrophysical flows - 
spherically symmetric or axisymmetric — transonic solutions pass through a critical 
point, it should be important to have a direct understanding of the nature of the 
critical points of the flow (and, of course, the physical solutions which pass through 
them). Without having to take recourse to the conventional approach of numerically 
integrating the governing non-linear flow equations, an alternative approach would be 
to adopt mathematical methods from the study of dynamical systems [19]. In accretion 
literature one would come across some relatively recent works which have indeed made 
use of the techniques of dynamical systems [H [20l [21] . Through this approach a complete 
and mathematically rigorous prescription can be made for the exclusive nature of the 
critical points in axisymmetric pseudo-Schwarzschild accretion, and, in consequence, 
the possible behaviour of the flow solutions passing through those points and their 
immediate neighbourhood. In conjunction with a knowledge of the boundary conditions 
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for physically feasible inflow solutions, this makes it possible to form an immediate 
qualitative notion of the solution topologies. What is more, this study has been shown 
to set a more restrictive condition on the specific angular momentum of sub-Keplerian 
flows passing through the critical points. 

While mentioning these general issues, it must be stressed however, that the specific 
purpose of the present work is to examine for the case of disc accretion some aspects of 
the feasibility of its transonic solutions, particularly their stability against the choice of 
a boundary condition and their long-time evolutionary properties on large length scales. 
To do so it would be necessary first to have a clear notion of at least the qualitative 
features of the phase portrait of the steady flow solutions, and its critical points. For the 
pedagogically simple and particular case of an axisymmetric flow driven by the classical 
Newtonian potential, it has been shown here that there are only two critical points for 
the steady flow, of which, for realistic boundary conditions, the outer one is a saddle 
point, while the inner one is a centre-type point. As an aside, it should be interesting 
here to note that this apparently simplistic scenario is identically reproduced for an 
axisymmetric accretion flow, under a fully rigorous general relativistic formalism, for 
certain values of the flow parameters [T8] . That these qualitatively physical conclusions 
about the flow could be drawn without falling back on the usual practice of a numerical 
integration of the steady flow equations, amply demonstrates the simplicity, the elegance 
and the power of the dynamical systems approach that has been adopted to study the 
thin accretion disc. 

Having gained an understanding of the nature of the critical points in the phase 
portrait, the question that is then taken up is about the preference of the accreting 
system for any particular velocity profile in the stationary limit of the flow, and a 
selection criterion thereof. In this regard the transonic solution has always been the 
favoured candidate. However, it is common knowledge from the study of dynamical 
systems, that a flow solution passing through a saddle point (the transonic flow in this 
case) cannot be realised physically [19]. To address this issue satisfactorily it must 
be appreciated that the real physical flow is not static in nature, but has an explicit 
time-dependence. With respect to this point it is tempting to subject the steady flow 
solutions to small perturbations in real time, and then study their behaviour. This has 
been done elsewhere [211 122] for the inviscid disc, and it has been shown that the steady 
inflow solutions of abiding interest are all stable under the influence of a linearised 
time-dependent perturbation on the mass inflow rate. 

Since one way or the other, no direct conclusion could be drawn about the selection 
of a particular solution through a perturbative technique, one could then try to have an 
understanding of a true selection mechanism and the attendant choice of a particular 
solution, by studying the evolution of the accreting system through real time. A model 
analog shows that it is indeed possible in principle for the temporal evolution to allow 
for the selection of an inflow solution that passes through the saddle point, and under 
restricted conditions it has been demonstrated that the selection criterion conforms to 
Bondi's minimum energy argument, which is invoked to favour the transonic solution 
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in the spherically symmetric case. Interestingly enough in this context, it has also been 
shown that although the perturbative study has offered no direct clue about the selection 
of a solution, the equation governing the propagation of an acoustic disturbance in the 
flow bears a close resemblance to the effective metric of an acoustic black hole. Use of 
this similarity has been made to argue that the flow would cross the acoustic horizon 
transonically. 



2. The equations of the flow and its fixed points 

It is a standard practice to consider a thin, rotating, axisymmetric, inviscid steady 
flow, with the condition of hydrostatic equilibrium imposed along the transverse 
direction [23l 124"] . The two equations which determine the drift in the radial direction 
are Euler's equation, 

v— + -— + <p'(r)-- = (1) 
dr p dr r A 

and the equation of continuity, 

4- (pvrH) = (2) 
dr 

in which <f>(r) is the generalised pseudo-Newtonian potential driving the flow (with the 
prime denoting a spatial derivative), A is the conserved angular momentum of the flow, 
P is the pressure of the flowing gas, and H = H(r) is the local thickness of the disc [24J, 
respectively. 

The pressure, P, is prescribed by an equation of state for the flow. As a general 
polytropic it is given as P = Kp 1 , in which K is a measure of the entropy in the flow 
and 7 is the polytropic exponent. The function H in ([2]) will be determined according 
to the way P has been prescribed [24], while transonicity in the flow will be measured 
by scaling the bulk velocity of the flow with the help of the local speed of sound, given 
as c s = {dP/dp) 1 ' 2 . 

With the polytropic relation thus specified for P, it is a straightforward exercise to 
set down in terms of the speed of sound, c s , a first integral of (pQ) as, 

v 2 \ 2 

^ + nc s 2 + 0(r) + — = £ (3) 

in which n = (7 — and the integration constant £ is the Bernoulli constant. Before 
moving on to find the first integral of (j2J) it should be important to derive the functional 
form of H. Assumption of hydrostatic equilibrium in the vertical direction deliver this 
form to be 

with the help of which, the first integral of (j^J) could be recast as 

cT n+1)V ~^ = t~ 2 M 2 (5) 

s 6' 47T 2 
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where Ai = { r yK) n fn [1] with rh, an integration constant itself, being physically the 
matter flow rate. 

To obtain the critical points of the flow, it should be necessary first to differentiate 
both (j3J) and ((SJ), and then, on combining the two resulting expressions, to arrive at 



d 2v 2 

(, 2 - pel) f(, 2 ) = *L 



A 2 xl (3 2 c 2 s 



3 — r — 

4>' 



(6) 



2 

with (3 2 = 2(7 + l) -1 . The critical points of the flow will be given by the condition 
that the entire right hand side of ([6]) will vanish along with the coefficient of d(f 2 )/dr 
on the left hand side. Explicitly written down, following some rearrangement of terms, 
this will give the two critical point conditions as, 



2„2 

sc 



r c <P'(r c ) - 



A 2 



3 — r, 



cf/(r c 



-1 -1 



(7) 



with the subscript c labelling critical point values. 

To fix the critical point coordinates, v c and r c , in terms of the system constants, 
one would have to make use of the conditions given by (j7]) along with (j3J), to obtain 



27 
7-1 



r c <P'{r c ) - 



A 2 



3 — r, 



<P"{r c ) 



+ 4>(r c ) + 



X 2 



2r 2 



£ 



(8) 



0'Oc) J 

from which it is easy to see that solutions of r c may be obtained in terms of A and £ 
only, i.e. r c = /i(A, £). Alternatively, r c could be fixed in terms of A and M.. By making 
use of the critical point conditions in (jHJ) one could write 

2(n+l) 



r c (j>'{r c ) - 



A 



3 — r. 



(9) 



with the obvious implication being that the dependence of r c will be given as r c = 
/2(A, .M). Comparing these two alternative means of fixing r c , the next logical step 
would be to say that for the fixed points, and for the solutions passing through them, 
it should suffice to specify either £ or M. pQ. 



3. Nature of the fixed points : A dynamical systems study 

The equations governing the flow in an accreting system are in general first-order 
non-linear differential equations. There is no standard prescription for a rigorous 
mathematical analysis of these equations. Therefore, for any understanding of the 
behaviour of the flow solutions, a numerical integration is in most cases the only recourse. 
On the other hand, an alternative approach could be made to this question, if the 
governing equations are set up to form a standard first-order autonomous dynamical 
system [19J. This is a very usual practice in general fluid dynamical studies [25], and 
short of carrying out any numerical integration, this approach allows for gaining physical 
insight into the behaviour of the flows to a surprising extent. As a first step towards 
this end, for the stationary polytropic flow, as given by (jSJ), it should be necessary to 
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parametrise this equation and set up a coupled autonomous first-order dynamical system 
as EEI 



■^(v 2 ) = 2v 2 



d 

d~T 

dr 
dr" 



A2 

— - rep H — 



3 — r — 



(10) 



in which r is an arbitrary mathematical parameter. With respect to accretion studies 
in particular, this kind of parametrisation has been reported before [U [20J, [21], but the 
present treatment, as far as disc accretion is concerned, will much more thoroughly 
highlight some serious questions about the feasibility of transonic flows. 

The critical points have themselves been fixed in terms of the flow constants. About 
these fixed point values, upon using a perturbation prescription of the kind v 2 = v 2 +Sv 2 , 
c s = c sc + ^ c s an d r = r c + 5r, one could derive a set of two autonomous first-order linear 
differential equations in the Sv 2 — 5r plane, with Sc 2 itself having to be first expressed 
in terms of 5r and 5v 2 , with the help of ([5]) — the continuity equation — as 

K _ 7-1 
4 " 7 + 1 

The resulting coupled set of linear equations in 6r and 5v 2 will be given as 



f5v 2 




A"(r c )" 


5r\ 




3 - r c - 










j)'(r c ) _ 





2v~l d7 



(Sv 2 



7 



+ 



A 

2 V7 + 1 

(3 2 clJ\r> 



5v 2 - 



2A 2 



r 



+ <t>'{r c ) +r c <P"(r c ) 



1 d \ 
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in which 



r c dr 



A 



cj)'(r c 



7 + 1 



-3, 



-Sv 2 



-B + 
A 



P 2 c 2 sc h-l 
7 + 1 



2 

7-1 
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A? 

r r 



dr 



B 



v 2 
r r 



- r, 



(12) 



0"(r c ) 



Trying solutions of the kind Sv 2 



r(r c ) <j)>{r c ) ' 
exp(fir) and 5r ~ exp(fir) in (fT2]) . will deliver the 



eigenvalues Q 



n- 



growth rates of 5v 2 and 5r 
2r c <P\r c )(3 2 cl 



as 



(7 + 1) 



(7 - 1) A - 2 7 (4 + A) + 2 7 £ ( 1 + - 



A- 



4 T +(7-l)i + 2 T Bh+y ^ 



(13) 



where Ak(^) is the local Keplerian angular momentum, expressed as A^(r) = r 3 0'(r). 

Once the position of a critical point, r c , has been ascertained, it is then a 
straightforward task to find the nature of that critical point by using r c in ( TL3l) . Since 
it has been discussed in Section [2] that r c is a function of A and £ (or M) for polytropic 
flows, it effectively implies that Q 2 can, in principle, be rendered as a function of the 
flow parameters. A generic conclusion that can be drawn about the critical points from 
the form of Q 2 in (fTBl . is that for a conserved pseudo-Schwarzschild axisymmetric flow 
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driven by any potential, the only admissible critical points will be saddle points and 
centre- type points. For a saddle point, Q 2 > 0, while for a centre- type point, Q 2 < 0. 
Once the behaviour of all the physically relevant critical points has been understood 
in this way, a complete qualitative picture of the flow solutions passing through these 
points (if they are saddle points), or in the neighbourhood of these points (if they are 
centre- type points), can be constructed, along with an impression of the direction that 
these solutions can have in the phase portrait of the flow [19]. 

A further interesting point that can be appreciated from the derived form of Q 2 , 
is related to the admissible range of values for a sub-Keplerian flow passing through a 
saddle point. It is self-evident that for this kind of flow, (A/Ak) 2 < 1 [6J. However, a 
look at ffl3|) will reveal that a more restrictive upper bound on A/Ak can be imposed 
under the requirement that Q 2 > for a saddle point, and this restriction will 
naturally be applicable to solutions which pass through such a point. This is entirely a 
physical conclusion, and yet its establishing has been achieved through a mathematical 
parametrisation of a dynamical system. 

To appreciate the practical usefulness of the method developed so far, it should be 
worthwhile to use the simple example of the Newtonian potential, = -GM/r as a 
special case. From (jSj) it will then be possible to find the spatial coordinates of the fixed 
points, which will be given by the two roots of a quadratic equation, whose solution will 
be given by 



(5 - 3 7 )GM 

10(7 - 1)£ 



1 ± 



10(5- 7 )( 7 -l)A 2 £ 



(14) 



(5-3 7 ) 2 (GM)2 

from which a conclusion that can be drawn is that for critical conditions to be obtained, 
there will have to be real solutions for r c , and this could only be achieved if the condition 

- c - < 5 - 3 ^ GM (15) 
Vl0(5- 7 )(7-l)f 

were to be maintained. 

The properties of these two critical points can be further analysed with the help of 
the eigenvalues of the stability matrix associated with the critical points. Upon using 
the Newtonian potential, <fi = —GM/r, in (fl3l) . it would be easy to arrive at 



n 2 = 2(3 2 c- 



7 \ GM 



7+1 



c A 



(16) 



A|(r c 

where C = (5 - 3 7 )(5 - 7 ) -1 and A|(r c ) = GMr c . 

The term in the square brackets in ffTB"]) will determine the sign of Q 2 . Knowing 
that two adjacent fixed points cannot be of the same nature [19] , i.e. they will differ 
in their respective signs of Q 2 , it is now evident that for the signs of all other factors 
in fflBI) remaining always unchanged, if r c is the length coordinate of the outer fixed 
point, for which (A/Ak) 2 < C> then Q 2 will be positive and the fixed point will be a 
saddle point, while if r c gives the inner fixed point, for which (A/Ak) 2 > C> then O 2 
will be negative and the fixed point will be centre-type. It is also remarkable that this 
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apparently simplistic scenario of two fixed points, with the outer one being a saddle 
point and the inner one being centre-type, can be identically reproduced in the case of a 
proper general relativistic flow on to a Schwarzschild black hole, under certain ranges of 
the flow parameters [18] . All the flows of physical interest in these cases are governed by 
the boundary condition that at large distances the drift velocity becomes vanishingly 
small, while the speed of sound approaches a constant value. This knowledge of the 
boundary condition, in conjunction with the nature of the fixed points, makes it possible 
to understand what the phase trajectories would look like. In a manner of speaking, 
the outer fixed point, which is a saddle point, may be dubbed the "sonic point" because 
one of the two trajectories passing through it is a transonic inflow solution, rising from 
subsonic values far away from the saddle point to attain supersonic values on length 
scales which are less than that of the saddle point. Transonicity is to be attained when 
the drift velocity equals the speed of the propagation of an acoustic disturbance, which, 
for this system is given by v2(7+ l) _1 ^ 2 c s [22]. Dwelling on this last point in somewhat 
greater detail, it may be noted that in some earlier works [21 E3] the speed of acoustic 
propagation has been scaled as an "effective" speed of sound. This scaling has been 
done by a constant factor that arises due to the chosen geometry of the thin disc flow 
(especially for a disc in vertical hydrostatic equilibrium, with the local disc height being 
a function of the acoustic velocity). This leads to the Mach number of the flow being 
scaled accordingly, and the critical condition for the flow will consequently be achieved 
when this scaled Mach number becomes unity. In the present treatment, on the other 
hand, the conventional definition of the Mach number has been adhered to, and as a 
result criticality will occur when the Mach number assumes the value \/2(7 + l) -1 / 2 . 
Either approach is entirely equivalent to the other, if one is mindful of the fact that the 
exact sonic condition (where the bulk flow velocity matches the local unsealed speed of 
sound) and the critical condition differ by a constant scaling factor only. This distinction 
disappears for the case of a spherically symmetric flow and for an axisymmetric flow 
with a constant disc height. On the other hand, the situation becomes radically more 
complicated for the case of a fully general relativistic flow. Here the critical Mach number 
is not a global constant, but is a local function of the critical point coordinates [18J, 
all of which, of course, makes it very much difficult a mathematical exercise to fix the 
critical point coordinates, calculate the slope of solutions passing through the critical 
points (in terms of physical relevance these points should be saddle points), and then 
numerically integrate these solutions under appropriate boundary conditions. 

It has been mentioned earlier that the particular properties of a saddle point in the 
phase portrait of thin disc flow solutions will also have a bearing on the possible range 
of values that the constant specific angular momentum, A, may be allowed to have. For 
solutions passing through the saddle point, it can be easily recognised from (jTj) that 
the flow will be sub-Keplerian [HE]. In that case the condition, (A/Ak) 2 < 1, shall 
hold good. In addition to this, for the specific case represented by ffT6l) . the saddle-type 
behaviour of the outer fixed point would also imply that there would be another upper 
bound on A, given by (A/Ak) 2 < C F° r the admissible range of the poly tropic index 7, 
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the possible range of £ would be < ( < 1/2. This would then imply that the latter 
bound on A would be more restrictive, as compared to the former. 

Saddle points are, however, inherently unstable, and to make a solution pass 
through such a point, after starting from an outer boundary condition, will entail an 
infinitely precise fine-tuning of that boundary condition [I] . This can be demonstrated 
through simple arguments. Going back to (fl2]) . the coupled set of linear differential 
equations in Sv 2 and Sr can be set down as 



in which the constant coefficients Qi, Q2, Q3 and Q4 are to be determined simply by 
inspection of (fT2"l) . It is also to be easily seen that Q\ = — <2 4 . This makes the integration 
of ( fT7l) a staightforward exercise that yields 



with C being an integration constant. Generally speaking ffl8|) is the equation of a conic 
section in the 8v 2 — 5r plane. If the origin of this plane were to be considered to have 
been shifted to the saddle point, then the condition for solutions passing through the 
origin, i.e. 5v 2 = 5r = 0, would be C = 0, which reduces (fT8|) to a pair of straight lines 
intersecting each other through the origin itself. All other solutions in the vicinity of 
the origin will, therefore, be hyperbolic in nature, a fact that is given by the condition 
(Q1 + Q2Q3) > 0. For the case of = —GM/r, this contention can be verified completely 
analytically, and this shows that even a very minute deviation from a precise boundary 
condition for transonicity (i.e. a boundary condition that will generate solutions to pass 
only through the origin, 5v 2 = 5r = 0) will take the stationary solution far away from a 
transonic state. This extreme sensitivity of transonic solutions to boundary conditions 
is entirely in keeping with the nature of saddle points. It may be imagined that in a 
proper astrophysical system such precise fulfillment of a boundary condition will make 
the transonic solution well-nigh physically non-realisable. Indeed, this difficulty, for any 
kind of accreting system, is readily appreciated by anyone trying to carry out a numerical 
integration of ([1]) to generate the transonic solutions, which can only be obtained when 
the numerics is first biased in favour of transonicity by using the saddle point condition 
itself as the boundary condition for numerical integration. 

There is another aspect to the physical non-realisability of transonic solutions, 
although this is somewhat mathematical in nature. Using the condition C = will 
make it easy to express 6v 2 in terms of 5r and vice versa. Going back to the set of linear 
equations given by ( TT2T) and choosing the second one of the two equations (the choice 
of the first would also have led to the same result), one gets 



which can be integrated for both the roots from an arbitrary initial value of Sr = [Sr]- m 
lying anywhere on the transonic solution, to a point 5r = e, with e being very close to 



d(5v 2 ) _ d (5v 2 ) /At 
d(6r) ~ d (Sr) /dr 



QiSv 2 + Q 2 5r 
Q 3 5v 2 + Q 4 5r 



(17) 



0,2 (Sr) 2 + 2Qi (5v 2 ) {5r) - Q 3 (5v 2 ) 2 + C = 



(18) 




(19) 
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the critical point given by Sr 

1 



0. With this it can be shown that 
d (5r) , 1 



Q2Q3 J[Sr] ln Sr 
from which it is easy to see that |r| - 



111 



(20) 



Q2Q3 

00 for e — > 0. This implies that the critical 
point may be reached along either of the separatrices, only after |r| has become infinitely 
large. This divergence of the parameter r indicates that in the stationary regime, 
solutions passing through the saddle point are not actual solutions, but separatrices of 
various classes of solutions [19J. This fact, coupled with the sensitivity of the stationary 
transonic solutions to the choice of an outer boundary condition, makes their feasibility 
a seriously questionable matter. 



4. Dynamic evolution as a selection mechanism for transonicity 

In Section [3] it has been demonstrated that generating a stationary solution through a 
saddle point (a transonic solution) will be impossible, physically speaking. Nevertheless, 
in accretion studies transonicity is not a matter of doubt. The key to this paradox 
lies in considering explicit time-dependence in the flow. This happens because the two 
stationary equations ([1]) and ([2]), as it is very much evident from their form, are invariant 
under the transformation v — ► —v (i.e. inflows and outflows are twin solutions of the 
same set of stationary equations), and it is this invariance that gives birth to the saddle 
point, which is an intersection point of the transonic inflow and outflow solutions. This 
invariance breaks down as soon as time-dependent terms are introduced in the governing 
equations, something that is very easy to see from 

dv dv 1 dP t A 2 
dt dr p dr r 3 
This will obviously imply that a choice of inflows (v < 0) or outflows (v > 0) has to be 
made right at the very beginning (at t = 0), and solutions generated thereafter will be 
free of all adverse conditions that one may associate with the presence of a saddle point 
in the stationary flow. 

In trying to make a time-dependent study, one could also assure oneself that a 
linearised perturbative analysis in real time cannot give any conclusive insight about 
the accreting system showing any preference for any particular solution. For inviscid, 
axisymmetric accretion, this fact has been clearly established [211 122] . One may then 
say that any selection mechanism based on explicit time-dependence has to be non- 
perturbative and evolutionary in character. 

However, even for the simple inviscid disc system, with explicit time-dependence 
taken into consideration, the equations of a compressible flow cannot be integrated 
analytically. Hence, to have any appreciation of the time-evolutionary selection of the 
critical solution, it should be instructive to consider an analogous model situation first. 
The model system being introduced here, describes the dynamics of a field y(x, t) that is 
analogous to a steadily accreting system with two fixed points, as has been discussed so 
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far. The details of the dynamic selection of the critical solutions in the model problem 
have been presented in the Appendix. 

To the extent that this model has been a good representative of the true physical 
situation, the whole argument for a time-dependent and non-perturbative method of 
selection developed in the Appendix, may now be extended to the actual problem of 
thin disc accretion. For a steady accretion disc many previous works have upheld the 
case for transonicity, although without explicitly addressing the issue of what special 
physical criterion may select the transonic solution to the exclusion of all other possible 
solutions. For the case of disc accretion on to black holes, Liang and Thompson [26] 
make a clear point by saying that "the solution for the radial drift velocity of thin disk 
accretion onto black holes must be transonic, and is analogous to the critical solution 
in spherical Bondi accretion, except for the presence of angular momentum." For the 
inviscid and thin disc at least, this is indeed a most crucial analogy, which will make it 
possible to invoke all the physical arguments used to uphold transonicity in spherically 
symmetric flows. 

Transonicity is a settled fact [3j [5] in spherically symmetric accretion, and this has 
been so because at every spatial point in the flow, velocity evolves at a much greater rate 
through time, in comparison with density, and this, therefore, excludes the possibility 
of matter accumulating in regions close to the surface of the accretor. As a result 
gravity wins over pressure at small distances and the system is naturally driven towards 
selecting the transonic solution. This sort of a situation will be more so obtained for the 
case of matter accreting on to a black hole, and it then raises the question of whether or 
not a likewise time-evolutionary mechanism should be at work for the selection of the 
critical inflow solution in an inviscid and thin accretion disc. This connection need not 
be so obvious, considering the fact that thin disc accretion involves, apart from a very 
different flow geometry, a whole range of different physical phenomena, not invoked for 
spherically symmetric flows. 

To have any analytical appreciation of how the temporal evolution acts as a selection 
mechanism in disc accretion, it should be instructive to prescribe the Newtonian 
potential for <p in (j2ip . as well as recast the pressure term in it by the polytropic 
relation, P = Kp 1 . All of this will lead to 



for which it should be noted that the question of generating transonicity in the inflow 
solution is not going to be affected overmuch by the specific choice of the simple 
Newtonian potential for 0, especially so in the vicinity of the far-off outer critical (saddle) 
point. It would be easy to see that to have a solution pass through this critical point, 
the temporal evolution should proceed in such a manner, that the inflow velocity would 
increase much faster in time than density at distances on the scale of the saddle point, 
where, for sub-Keplerian flows, gravity, going as r~ 2 , dominates the rotational effects, 
which go as r~ 3 . 








(22) 
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But what should be the key physical criterion guiding this dynamic selection of the 
transonic solution? It is to be stressed once again that the selection principle would very 
likely be the same as in the spherically symmetric case, where the transonic solution is 
chosen by dint of its corresponding to a configuration of the lowest possible energy [5]. 
To have an understanding of this, an analytical treatment of ( 1221) may be carried out, 
with the density gradient being neglected as a working approximation. On large length 
scales, where the flow is highly subsonic, this is an especially effective approximation, and 
it will allow for treating the evolution of the velocity field to be largely independent of 
the density evolution. The resulting non-linear partial differential equation for velocity 
is then to be integrated by the method of characteristics [27] (the mathematical details 
of this method have been outlined for the model problem presented in the Appendix) 
to get a solution that can be written as 

v 2 GM A 2 



2 r 2r 2 



( vr at 

exp 



(23) 



1 + v/aj \ar r 0/ 

in which r = GM/a 2 , and a itself is an integration constant deriving from the spatial 
part of the characteristic equations, and, therefore, should, in general, be a function of 
time. The form of the function T is to be determined from the physically realistic initial 
condition v — at t — for all r, which will render T as 

GMi \ 2 i 2 

*° = + io^o 1 " (24) 

On examining the argument of T in ( 1231) and then studying its long-time behaviour, it 
will be seen from ( 12^1) that for t — > oo, the selected solution will correspond to the 
condition 

v 2 GM X 2 



2 r + 2^=°- < 25 > 
It can now be seen that prior to the evolution, the system had no bulk motion and 
that the radial drift velocity was given flatly everywhere by v = 0. This, of course, gives 
the condition that initially the total specific mechanical energy of the system was zero. 
Then at t = both a gravitational mechanism is activated in this system and some 
angular momentum is imparted to it. This will induce a potential —GM/r everywhere, 
and at the same time start a rotational motion, respectively. The system will then start 
evolving in time, with the velocity v at each point in space evolving temporally according 
to (T231 . Finally the system will restore itself to a steady state in such a manner that the 
total specific mechanical energy at the end of the evolution (for t — > oo) will remain 
the same as at the beginning (at t — 0), a condition that is given by (|25|) . whose left 
hand side gives the sum of the specific kinetic energy, the specific gravitational potential 
energy and the specific rotational energy, respectively. This sum is zero, and therefore, 
under the given initial condition, this must be the steady state corresponding to the 
minimum possible total specific energy of the system. Hence, this is the configuration 
that is dynamically and non-perturbatively selected. It is now conceivable that if the 
pressure term were to be taken into account, then the solution that would be dynamically 



Evolution of transonicity in an accretion disc 



13 



selected, would be the one that would pass through the saddle point, since, as Bondi 
had analogously conjectured for spherically symmetric accretion [3], this would be the 
one to satisfy the criterion of minimum energy. 

Interestingly enough, the issue of the primacy of transonicity can also be addressed 
from a very different, and somewhat unlikely, perspective. It has been mentioned earlier 
that carrying out a linearised perturbative analysis in real time on stationary flows will 
not indicate any particular solution to be favoured over all the others. This line of 
thinking may now be subjected to a closer scrutiny. 

Under the condition of hydrostatic equilibrium in the vertical direction, the time- 
dependent generalisation of the governing equations for an axisymmetric pseudo- 
Schwarzschild disc is given by (T2~T1) . as well as 



d£ 1 d 
dt r dr 



Evr) = 



(26) 



in which the surface density of the thin disc S, is to be expressed as S = pH 
Making use of (jl]) and the polytropic relation P = Kp 1 , (T2"6"l) can be rendered as 



ii 3 1 '2 (^ , y 



p 



(7+l)/2 



,3/2 



V- 



0. 



(27) 



Defining a new variable / = p^ 7+1 ^ 2 t>r 3 / 2 /v / 7 , it is quite obvious from the form of 
( f2"T|) that the stationary value of / will be a constant, fo, which can be closely identified 
with the matter flux rate. This follows a similar approach to spherically symmetric 
flows established in earlier works [281 [29]. The present treatment, of course, pertains to 
a disc flow being driven by a general pseudo- Newtonian potential, <f)(r). In this system, a 
perturbation prescription of the form v{r, t) = vo(r)+v'(r, t) and p(r, t) = po(r)+p'(r, t), 
will give, on linearising in the primed quantities, 



L 

fo 



7 + 1 



Po 



V' 
v 0. 



(28) 



with the subscript denoting stationary values in all cases. From (1271) . it then becomes 
possible to set down the density fluctuations p', in terms of /' as 

'df r 



dp' n2 V OPO 

dt +P fo 



dr 







(29) 



as before. Combining (1251) and (12"9"|) will then render the velocity 



with P 2 = 2(7 + 1) 
fluctuations as 

dv' vq 
~dt = To 

which, upon a further partial differentiation with respect to time, will give 



dt dr 



<9V 



d_ 

dt 



fo 



df 
dt 



d 

+ dt 



fo 



df 
dr 



From (12T|) the linearised fluctuating part could be extracted as 
dv' 



dt 



d_ 

dr 



I 2 P 
V V +C s0 — 

Po 







(30) 



(31) 



(32) 
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with c s0 being the speed of sound in the steady state. Differentiating (132]) partially with 
respect to t, and making use of (1291) . fl30l and fl3TT) to substitute for all the first and 
second-order derivatives of v' and p', will deliver the result 



0_ 

dt 



fo 



df 

dt 



+ 







Ot 



fo 



df 
dr 



+ 







+ 



dr 
d_ 
dr 



fo 
fo 



df 

dt 



df 

dr 



(33) 



all of whose terms can be ultimately rendered into a compact formulation that looks 
like 

d tl (^d u f) = (34) 

in which the Greek indices are made to run from to 1, with the identification that 
stands for t, and 1 stands for r. An inspection of the terms on the left hand side of (J33J) 
will then allow for constructing the symmetric matrix 

f ^ = ( 1 

fo \vo 

Now in Lorentzian geometry the d'Alembertian for a scalar in curved space is given 
in terms of the metric g Mi , by [301 

1 



v o — (3 c s0 



(35) 



(36) 



with g^ u being the inverse of the matrix implied by g M!/ . Comparing fl34l) and fl36l) it 
would be tempting to look for an exact equivalence between and \f—gg^ v . This, 
however, cannot be done in a general sense. What can be appreciated, nevertheless, is 
that gives an equation for /' which is of the type given by ( 1361) . The metrical part 
of (l34"j) . as given by fl35l) . may then be extracted, and its inverse will incorporate the 
notion of a sonic horizon of an acoustic black hole when Vq = f3 2 c^ Q . This point of view 
does not make for a perfect acoustic analogue model, but it has some similar features to 
the metric of a wave equation for a scalar field in curved space-time, obtained through 
a somewhat different approach, in which the velocity of an an irrotational, inviscid 
and barotropic fluid flow is first represented as the gradient of a scalar function ip, i.e. 
v = — V?/>, and then a perturbation is imposed on this scalar function (3TJJ, [3Tj . 

The discussion above indicates that the physics of supersonic acoustic flows closely 
corresponds to many features of black hole physics. This closeness of form is very 
intriguing, as well as instructive. For a black hole, infalling matter crosses the event 
horizon maximally, i.e. at the greatest possible speed. By analogy the same thing may 
be said of matter crossing the sonic horizon of a fluid flow. Indeed, this has been a 
long-standing conjecture for the case of spherically symmetric accretion on to a point 
sink [3], [5]. That this fact can be appreciated for the accretion problem through a 
perturbative result, as given by fl33|) . is quite remarkable. This is because conventional 
wisdom would have it that one would be quite unable to have any understanding of the 
special status of any inflow solution solely through a perturbative technique [5]. It is 
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the transonic solution that crosses the sonic horizon at the greatest possible rate, and 
the near similarity of form between (133]) and fl36l) may very well be indicative of the 
primacy of the transonic solution. If such an insight were truly to be had with the help 
of the perturbation equation, then the perturbative linear stability analysis might not 
have been carried out in vain after all. 

5. Concluding remarks 

The difficulties against transonicity, arising from the choice of purely stationary 
equations, have been addressed through non-perturbative dynamics. This has been 
a relatively simple exercise to carry out in the Newtonian construct of space and time. 
It may be readily appreciated, however, that a rigorously general relativistic flow will 
not lend itself so easily to such a treatment as has been executed here for a pseudo- 
Newtonian system. Nevertheless, it should be helpful to note that particular pseudo- 
Newtonian systems [32] bear a very close resemblance to an actual general relativistic 
system [Hj. This closeness is a reason to believe that the non-perturbative approach to 
the selection of transonicity demonstrated in this paper may also be very relevant for 
general relativistic flows. 

As a final point, it may be mentioned that realistically speaking, a disc system 
should involve viscosity as a means of transporting angular momentum, to make infall 
a sustained global process. However, in that event, viscosity will also bring about 
dissipation of energy in the flow, something that will make the minimum energy criterion 
an ineffective instrument for identifying the transonic solution. Besides this, viscosity 
will violate Lorentzian invariance, which is very crucial for constructing an analogue 
gravity model. With the breakdown of the Lorentzian invariance, the quest for an 
acoustic analogue of a black hole becomes a most difficult one, especially when the fluid 
system under study is compressible in nature. For incompressible flows, on the other 
hand, some studies have considered the issue of an analogue gravity model vis-a-vis 
viscosity in a semi- quantitative manner. In analysing the outflow of a very shallow layer 
of water on a flat surface, it has been discussed that the formation of an abrupt hydraulic 
jump in the flow is entirely due to viscosity, and many features of the jump itself can 
be closely connected to an acoustic white hole [331 El] • But more importantly the basic 
properties of surface gravity waves in a shallow layer of water remain unchanged in the 
vicinity of the jump, inspite of viscosity [35J. Going by this analogy, it is conceivable 
that the properties of an acoustic wave in a compressible flow (such as the accretion 
process is) will likewise remain unaffected. This argument, however, will not make it 
directly possible to establish precise relations like (1341) and ( l35l) . which give a clear 
mathematical result to indicate that transonicity might be favoured. Yet transonicity 
will very likely be achieved because regardless of energy dissipation due to viscosity, the 
stationary matter flux rate will be conserved and the flow will continue to proceed at 
the greatest possible rate. 
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Appendix 

The model system that has been introduced here serves to show that non-realisable 
separatrices passing through the saddle point in the stationary regime, can behave like 
proper physical flows, when the dynamics, as opposed to the statics, is to be followed. 
The dynamics of the field y(x, t) is described as 

^ + {y ~ x) ^ = y + 2x + x2 {AA) 
whose static limit leads to 

dy y + 2x + x 2 

dx y — x 

and which, viewed as a dynamical system, is seen as 

d V ,0,2 

— = y + 2x + x 

dr 

dx , . . 

— =y-X. (A.3) 

In the y — x space, the fixed points (x c ,y c ) are to be found at (0,0) and (—3,-3). 
A linear stability analysis of the fixed points in r space gives the eigenvalues Q, by 
Q 2 = 1 + 2(x c + 1). It is then easy to see that (0, 0) is a saddle point while (—3, —3) 
is a centre-type point. This mathematical model is very much similar to the stationary 
inviscid accretion disc driven by the Newtonian potential, with two critical points, of 
which the outer one is a saddle point, while the inner one is centre-type. The integral 
curves of the model system are obtained from (1A.2I) as 

2 

y 2 - 2xy - 2x 2 x 3 = c (A.4) 

3 

with the solutions passing through the saddle point being given by c = 0. 

To explore the temporal dynamics, and to obtain a solution to flA.ll) . it would be 
necessary to apply the method of characteristics [27J. This involves writing 

^ = = (A.5) 

1 y — x y + 2x + x 1 

The subsequent task is to find two constants c\ and c<i from the above set and the 

general solution of (lA.ljl would then be given by c\ = F(c2), where the function F is to 

be determined from the initial conditions. It is easy to see that one of the constants of 

integration is clearly the c of flA.4j) . Hence, writing c\ = c, and using flA.4j) in the first 

part of ( 1A.5I) . will give 



dt = ± I dx (A.6) 

' v/3x 2 + (2/3) x 3 + c 
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which solves the problem in principle. To put this in a usable form, the integration in 
flA.6j) would have to be carried out. For small x (the most important region, since it is 
near the saddle), the x 3 term may be left out to a good approximation. Further, only 
the positive sign in the right hand side of (1A.6I) is to be chosen by the physical argument 
that the system is to evolve through a positive range of t (time) values. Integration of 
flA.6j) will then lead to the result 

-\/3t 



(x + \/x 2 + e~ w "" = c 2 
which will then make the solution of flA.lj) look like 



y 2 - 2xy - 2x 2 — -x 3 = F 



x + 



(y - x) 



-r 
9 



(A.7) 



(A.8) 



The evolution of the system is to be followed from the initial condition that y = 
for all x at t = 0. Dropping the x 3 term again for small x, will allow for determining the 
form of the function F as F(z) = —3(2 — \f?>)z 2 . The solution, consequently, becomes 

2 



y 2 -2xy- 



-2x 2 — x 3 
3 



(2-V5 




(y - %y 



—X" 

9 



,-2V3t 



(A.9) 



and as t — >oo, the steady solution that would be selected would be 



y 



2xy — 2x 2 



-x 
3 







(A.10) 



which is actually the equation for the separatrices. It is worth stressing the remarkable 
feature of this result. The evolution started under conditions far removed from 
transonicity. In fact, it started as y = for all x (in the vicinity of the origin of 
coordinates) at t — 0. The evolution proceeded through a myriad of possible steady 
state solutions (all arguably stable under a linear stability analysis) and then in the 
stationary limit, selected the separatrices. This is a convincing demonstration that it 
is in principle possible for apparently non-realisable separatrices in the steady regime, 
to become eminently realisable physically, when the temporal evolution of the system is 
followed. 
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